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METHOD OF ESTIMATING A SIGNAL-TO- 
INTERFERENCE+NOISE RATIO (SINR) 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to wireless communication, and more particularly, a 
method of estimating a signal-to-interference+noise ratio. 

2. Description of Related Art 

Signal-to-Interference+Noise Ratio (SINR) is an important metric of 
communication link quality. SINR estimation is of particular importance for wireless 
data systems where resources are shared dynamically amongst users. Some applications 
of SINR estimates are: a) Power Control in CDMA Systems: the receiver estimates the 
SINR, compares it to a target and commands the transmitter to increase/decrease its 
transmitted power; and b) Rate Adaptation: the information bit-rate assigned to a user 
can be dynamically varied based on its link quality and the system load. While such 
adaptation has limited use in voice systems, it is extremely useful for wireless data 
systems. Consequently, inaccurate SINR estimates can severely degrade performance 
and resource utilization. 

SUMMARY OF THE INVENTION 

In the method of estimating a signal-to-interference+noise ratio (SINR) according 
to the present invention, a conventionally known SINR estimate is scaled and translated 
to produce a SINR estimate that has a substantially reduced mean square error. The 
SINR estimator according to the present invention operates on either pilot symbol 
samples or data symbol samples. Further according to the present invention, a pilot 
symbol based SINR estimate and a data symbol based SINR estimate are weighted and 
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combined to generate a composite SINK estimate having a reduced mean square error as 
compared to either of the individual estimators. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 

The present invention will become more fully understood from the detailed 
description given herein below and the accompanying drawings, which are given by way 
of illustration only, and thus are not limitative of the present invention, and wherein: 

Figure 1 illustrates pilot and data symbols multiplexed in a time slot; 
10 Figure 2 illustrates the gain from the SDSTR estimator according to an embodiment 

of the present invention as compared to a conventional estimate based on pilot symbols 
only for different smoothing factors. 

Figure 3 illustrates optimal weights assigned to a pilot-symbol based estimator 
and a data-symbol based estimator in a combined estimators; and 
15 Figure 4 illustrates the performance of the combined estimator using the 

estimated weights. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

20 In describing the method of estimating the signal-to-inteference+noise ratio 

according to the present invention, only Binary Phase Shift Keying (BPSK) modulation 
is considered although the methods and related analysis can be extended to other 
signaling schemes. Noise and interference are modeled together as additive white 
Gaussian noise (AWGN), but as will be appreciated from the following disclosure, this 

25 should not limit the application of the method according to the present invention. 
Transmission is organized into fixed duration timeslots, each containing pilot and data 
symbols as shown in Fig. 1 . The channel attenuation and phase shift, assumed fixed over 
a timeslot, are treated as unknown constants rather than as random variables (no fading 
assumption). 

30 First, to better understand the method according to the present invention, the 

conventional SINR estimation will be described. Typically, the received signal 
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corresponding to the 7'th transmitted symbol (pilot or data) in the Ath time slot is defined 
as 

(1) 

hj=%Mk+Ski j=\a,...,N, 

5 

where /Uk represents the received signal amplitude (product of transmitted amplitude and 
channel gain), ^-^^ is a random variable that represents the noise+interference, 
represents the symbol-value, and N is the number of samples (pilot or data). Information 
symbols could be +1 or -1 (in BPSK), while it is assumed (without any loss of 
10 generality) that pilot symbols are always +1. It is also assumed that the distribution that 
characterizes the noise+interference is Gaussian with zero mean and variance cr^. The 
SINR in the kth time slot is defined as: 

(2) 

CT 

15 and is the parameter to be estimated. 

The groups of N sample points (data or pilot) could correspond to a time slot in 
CDMA systems or frames in TDMA systems. A well-known pilot-symbol sample based 
estimator of SINR is computed as the ratio of the square of the sample mean of the 
received to the pilot-symbol sample Y (based on N sample points in a group) to the 

20 corresponding sample variance is a reasonable estimator of the SINR. Estimators based 
on this ratio are called Squared Mean By Variance or SMV estimators. Different SMV 
estimators have been studied in the literature and they only differ in the normalization 
constant used to compute the sample variance. 

For the case where the {Yk,} values correspond to pilot symbols, define the 

25 sample mean and unbiased sample variance for the Ath time slot as 

(3) 

_ (4) 
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Then, 0^ ={Yt)^ /Sl is one possible SMV estimator (commonly used). The Maximum 
Likelihood (ML) estimator of the SINR is also an SMV estimator where a biased sample 
variance estimate is used in the ratio (normalization is by N rather than N-1). The signal 
power, jul , varies due to channel fading. However, the noise variance, a^, changes very 

5 slowly with time, typically with the addition (departure) of a call; therefore, one can 
improve the overall quality of the SINR estimate by using a longer-term moving average 
estimate of the sample variance. This makes the "effective" number of samples used in 
the sample variance estimate larger and therefore more accurate. One simple method to 
accomplish this is through exponential smoothing of a set of sample variance estimates 

10 (called Exponentially Weighted Moving Average or EWMA). The smoothed sample 
variance estimate through the Ath time slot would be given by: 



where r is the smoothing factor determined according to desired design parameters 
15 andO < r<\. The SINR estimate at the end of k time slots then becomes: 



The accuracy (mean and variance) of SMV estimators based on using a single group of 
20 pilots with N sample points is known in the art. However, aside from simulation results, 
the accuracy of SMV estimators that utilize EWMA for sample variance estimation is not 
known. 

The inventors have analytically analyzed the accuracy of SMV estimators that 
utilize EWMA for the sample variance estimate and have developed accurate 
25 approximations for the bias, variance and MSE of such estimators. From this analysis, 
the inventors developed a characterization of the bias and variance, and hence the MSE, 
of such estimators. The inventors, based on their analysis, developed a scaling and 
translation methodology for the SINR estimators that substantially reduces the MSE. The 
scaled and translated EWMA SINR estimator according to the present invention is: 



(5) 



al={\-r)6rl_^+rSl k>\. 



(6) 
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(7) 



where is the translation and n is the scahng factor defined as: 

5 (8) 

,^_ (A^-l)(2-r) _ 

Note that this scaling factor ?j depends only on parameters that are fixed and known at 
the receiver and do not require any separate estimation procedure. 

The gains due to scaling and translation (measured as the reduction in mean- 

10 squared error) in using 0^ instead of S/, can be substantial especially for larger values of 
the smoothing factor r as seen in Figure 2. Figure 2 illustrates the gain (measured as the 
reduction in MSE) based on pilot symbols only for different smoothing factors r , and 
show that even for values of r around 0.5-0.6, the gain is around 30-40%. 

SMV estimators that use data symbols are also known in the art. A commonly 

15 used SMV estimator based on data symbols, called the non-coherent estimator, is one 
where the sequence {Ykj} is replaced by the sequence of its absolute values i.e. 
{Zkj=\Ykj\}. The SnSTR estimate for the Ath time slot is the ratio of the sample mean and 
sample variance of the sequence {Zt,}. Smoothing of the sample variance of the {Zig} 
sequence via the EWMA approach may also be used to improve accuracy. While the 

20 scaling and translation approach developed in this invention can be extended to any 
SMV estimator based on data symbols, including the non-coherent estimator, the 

preferred approach is outlined in concurrently filed Application No. 

entitled "Method of Estimating a Signal-to-Interference+Noise Estimate Using Data 
Samples" by the inventors of the subject application. 

25 The accuracy of the SINR estimate can be further enhanced by simultaneously 

ufilizing estimates from pilot symbols as well as data symbols in the estimation process. 
As mentioned before pilot symbols have the advantage of being known at the receiver 
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but are relatively fewer in number. Data symbols are more plentiful than pilot symbols 
but they are unknown at the receiver. 

First, a method for the general problem of optimally linearly combining two 
estimates of an unknown parameter will be described. With that foundation, the 

5 combining method will be applied to SINR estimation with pilot and data symbols. 

Suppose there are two noisy estimates, ©, and ©2, respectively, of a parameter 
©. Furthermore let the biases of the two estimates be B^ = E[@,]-@ and B2 = E[&J-@ , 
respectively, and the variances be rf and r^, respectively. The MSB of 0, is 
MSE, ==Bf+Tf, i = l,2. Consider the class of estimators formed by an arbitrary linear 

10 combination of 0, and ©2 



(9) 

0 = ai0i + a2©2. 

Minimum mean squared error weights and ^2 are obtained by minimizing 
15 E(©-©)^. The following are the minimum mean squared error weights 

(10) 

ar= (^±A)^^-^© 

' (© + B^y+ rl (© + 5i )' + Tlr\ 



:f (0 + B^f+ x\ (0 + 5, + zIt\ 

(11) 



The corresponding minimum MSE is 

(12) 

MSE""' = r— ^ ^ ^ . 

1 + (0 + 5,)' ItI + (0 + ^2 ) 

Unfortunately, the optimum weights to use depend on the value of © . Later in 

this section, we shall describe how the weights can be estimated in practice. 



-6- 



Docket No.: 2925-555P 
Lucent Ref. No.: 122646/Jeske 6-11 



When two unbiased estimators are being combined (i.e. B]=0 and 82=0) and it is 
desired that the combined estimators also be unbiased, then the weights need to be 
constrained to sum to unity. In this case the optimal weight on ©, can be shown to be 

(13) 

MSE2 - B^B^ 

^ " ~ MSE^ + MSE^ - IB^B^ 

and the optimal weight on ©2 is thus \-a = {,MSE^-BiB2)l{MSE^ + MSE^-IB-^B^). The 
MSB of the optimal constrained estimator is 

(14) 

. rr.^on, MSE, USE. - B} bI 
MSE ' = ' — - — . 

10 When both estimators are unbiased Equation (13) reduces to the known result for 
optimally combining (linearly) two independent unbiased estimates. 

The preceding results will now be applied to the problem of SINR estimation. 
Consider two estimators of the SfNR, . For the purpose of illustration we select the 
smoothed, scaled and translated pilot-sample based estimator, say ©f , and the smoothed, 

15 scaled and translated data-sample based estimator (e.g. the non-coherent estimator), say 
e>f , as the two estimators to be combined. The inventors have shown that the biases and 
variances of these estimators are approximated as follows 



5f 

..y,2(^r (l.^-.©.„)- ^i.2A^©„j _ 

20 

In Equation (15), the superscript and subscript P denotes the pilot sample based 
variables, the superscript and subscript D denotes the data sample based variables, No 
and Np are the number of data and pilot symbols respectively per time slot, 
={Np -l)(2-r)/r , 7^ = (^o -l)(2-r)/j- (where r is the EWMA smoothing factor) and the 
25 quantity 
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where 




(erfc in the above equation is the well-know complementary error function). Using 
Equation (15) with either Equation (10) and Equation (11) or Equation (13) gives the 
optimal combined estimators of and ©f. Since the biases and variances, and 
therefore the weights themselves depend on ©j, which is unknown, one can only 

10 estimate the weights in practice. An estimate of the weights can be obtained by replacing 
by ©f in the corresponding expressions. One could, alternatively, replace ©^ by ©f 
to obtain estimates of these quantities. However, since ©f is an unbiased estimate of , 
the accuracy of the composite estimate is better, when ©[ is used. Note that the above 
technique for combining may be used with any two estimators as long as the appropriate 

15 bias and variance expressions are used in the weight computation. For example, one can 
combine unsmoothed, scaled and translated pilot and data-sample based estimators using 
the same approach. 

Let CI denote the combined estimator of ©f and ©f that has constrained weights 
and C2 denote the combined estimator with unconstrained weights. The optimal weights 

20 used when ©f and ©f are combined for CI and C2 are shown in Figure 3, as a function 
of . In the figure-legend BCSV (BCSV-Z) stands for the smoothed, scaled and 
translated pilot (data) sample based estimator. The root mean squared error (RMSE) 
values for CI and C2 when estimated weights are used are shown in Figure 4. For 
comparison, the corresponding RMSE values for each of the component estimators, 

25 BCSV and BCSV-Z, are also shown. (The label "A" in the figure merely means that 
numerical results were obtained from analytic formulas rather than simulations). Despite 
the error in estimating the weights, both combined estimators outperform the individual 
estimators for the entire range of SINK considered. 
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The invention being thus described, it will be obvious that the same may be 
varied in many ways. Such variations are not to be regarded as a departure from the 
spirit and scope of the invention, and all such modifications as would be obvious to one 
skilled in the art are intended to be included within the scope of the following claims. 



